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ABSTRACT 

We use published spectroscopic and photometric data for 8 Type la supernovae to 
construct a dispersion spectrum for this class of object, showing their diversity over the 
wavelength range 3700 A to 7100 A. We find that the B and V bands are the spectral 
regions with the least dispersion, while the U band below 4100A is more diverse. Some 
spectral features such as the Si line at 6 150 A are also highly diverse. We then construct 
two objective measures of 'peculiarity' by (i) using the deviation of individual objects 
from the average SN la spectrum compared to the typical dispersion and (ii) applying 
principle component analysis. We demonstrate these methods on several SNe la that 
have previously been classified as peculiar. 
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. 1 INTRODUCTION 

Type la supernovae (SNe la) are observed to display re- 
markably homogeneous light curves in broadband photom- 
etry. They are currently classified according to a single- 
parameter fa mily of peak l uminosity as a function of light 
curve shape (iPhillii: which allows them to be used 

as accurate cosmological distance indicators. 

However, we know that this one-parameter family does 
not completely encompass the variety of SNe la observed. 
Once we have reduced any contribution from measurement 
error to the point where its effects are negligible the intrinsic 
diversity of SNe la themselves represent the limiting accu- 
racy to which a one-parameter family can aspire. 

This diversity manifests itself as differences in the SN 
la spectra. It is important to assess the diversity in partic- 
ular features of SN la spectra and find those hot spots that 
have the most significant effect on the broadband photome- 
try used for cosmology. Identifying these points of variabil- 
ity in SN la spectra has a multi-fold benefit. (1) It allows 
us to identify which spectral regions are the most stable for 
photometry and thus design future dedicated SN cosmology 
probes to utilise these stable wavelength regions. (2) It al- 
lows us to classify SN la based on certain spectral features 
into sub-types that have a tighter magnitude-redshift rela- 
tion than the entire SN la population. (3) It allows us to put 
weighted error bars on the points in the magnitude-redshift 
relation depending on which part of a redshifted SN la spec- 
trum was sampled by the filter set in use. 

E-mail: jbjames@pliysics.usyd.edu.au (JBJ); tama- 
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The quest to understand the diversity in the type la 
supernova population is accelerating, with many groups in- 
vestigating this in an attempt to control systematic un- 
certainties associated with using type la supernovae as 
standard candles. Empirical methods for studying the va- 
riety in the type la population focus on the linked as- 
pects of photometric diversity ((jha et al.l l2005l) and spec- 
troscopic diversity. Many of the spectral studies concen- 
trate on fitting spectral lines (such as the Si II feature), 
measuring how they evolve with time, relating these to 
physical parameters such as metallically or eject a veloc- 
ity, and studying correlations between these features and 
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peak magnitude or light curve width (iNuerent et al. 




iFiliDDenkcl Il997l iBenetti et al.l 120041 iBlondin et al. 
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"Gallagher et al."'2005; 


'Benetti ela]J|2005V Efforts to 



compare these observational parame ters to the results of 
theoretical models are ongoing (e.g . Ilravag lio et al.l l2005l : 
iBranch et alJl200lEoO^ : lBongard et al.l l2006h. and this im- 



portant topic promises to allow the prediction of potential 
changes in type la supernova demographics over time. In 
the meantime observational correlations, such as the width- 
luminosity relation, remain the strongest methods for con- 
trolling the effects of supernova diversity on cosmological 
parameter estimation. 

In this paper we design an empirical measure of diver- 
sity that uses a new approach. Instead of fitting spectral 
lines we subject the entire spectrum to a number of differ- 
ent statistical analyses. These show promise for providing an 
objective method of distinguishing peculiar supernovae (at 
the very least) and for finding correlations amongst subsets 
of the type la population. In Section |21 we describe the data 
set and our analysis technique. In Section El we provide plots 
of the spectral variability of SNe la as a function of wave- 
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Table 1. Catalogue of SNe la spectra used in this study. 



SN 


epoch (days) 


Amin (A) 


Amax (A) 


Reference' 


1981B 





1272 


8392 


Br83 


1989B 


— 1 


3001 


11167 


We94 


1992A 


— 1 


3564 


7100 


Ki93 


1994D 


—2 


3402 


9132 


Pa96 







3069 


10130 


Br 




2 


3466 


9204 


Pa96 




2 


4439 


7013 


Pa96 


1994S 





3120 


11300 


Br 


1996X 


-2 


3748 


9917 


SalOl 







3079 


10670 


SalOl 




1 


3078 


10669 


SalOl 


1998aq 





3720 


7500 


Br03 




1 


3720 


7511 


Br03 




2 


3720 


7521 


Br03 


1999ee 


-2 


3262 


9983 


Ha02 




-1 


3559 


9134 


Ho 







3262 


9687 


Ha02 



^References.— Br83: Branch et al. 1983; Br03: Branch et al. 
2003; Ha02: Hamuy et al. 2002; Ki93: Kirshner et al. 1993; Pa96: 
Patat et al. 1996; SalOl: Salvo 

length and discuss the features. In Section 0] we provide a 
prescription for determining the degree of 'peculiarity' of a 
supernova based on its spectrum near maximum light. 



2 DATA PREPARATION 

The sample consists of 18 previously published spectra taken 
for 8 low redshift type la SNe. Table Q lists the spectra and 
the source catalogues. These represent the subset of pub- 
lished objects that are confirmed to be normal SNe la and 
have good quality spectra within two days of maximum light 
as well as established Amis values, B magnitudes at maxi- 
mum light (Mb), extinction values {E{B — V)) and distance 
moduli (/io)- Table 121 lists the photometric properties of the 
source objects. 

Correction for dust ext inction is made according to the 
standard extinction law o feavage Mathis (1979), usmg 
colour excess values E{B — V) from IPhillips et al.l (| 199911 
where available. 

It is a familiar task to convert apparent magnitudes 
to absolute magnitudes in order to compare the brightness 
of objects at different distances. In order to compare the 
18 spectra in our sample we have to perform the spectral 
equivalent of converting from apparent to absolute magni- 
tude. We do this by multiplicatively scaling the spectra so 
they give the correct absolute magnitude in the B or V band 
after synthetic photometry. 

There is more than one way to determine the 'correct' 
absolute magnitude. We could choose to scale each spectrum 
to the broadband magnitude of its source supernova, once 
corrected for extinction, distance modulus and Ami5(B). 
However the observational uncertainty on nearby SN la mag- 
nitudes is typically about 0.2 magnitudes, while the stan- 
dard deviation in magnitude values is typic ally about 0.4 
and .2 magnitudes in B and V respectively (|Hamuv et al.l 
Il996l l. This makes it difficult with the current data to dis- 



tinguish broadband magnitude variations due to intrinsic 
spectral variability from variations due to observational er- 
ror. 

Instead we choose to scale each spectrum to a fixed 
broadband magnitude. This allows us to probe the spectral 
diversity without contamination from observational magni- 
tude uncertainties. When we scale to a fixed magnitude we 
artificially minimise the diversity in that region. To examine 
the magnitude of this effect, we scale the spectra to fixed 
Mb in the first instance and to fixed My in the second. 
The choice of scaling has a small effect compared to the 
size of the spectral variations. The mean values of the B 
and V absolute magnitudes are Mb = —19.16 ± 0.05 and 
Mv = -19.1 1 =b 0.02 (Icr) agreei ng: with the better perfor- 
mance of the lHamuv et all (^9^) relation in the V band. 

Scaling in a single band is achieved by allowing for a 
linear relationship between the synthetic photometric mag- 
nitude and the logarithm of the scaling factor (luminosity)^, 

M = a log /c + 6, (4) 

The value of the slope constant a is constant for all objects 
within a particular band. The offset b varies depending on 
the units of the spectrum, and is determined for each object 
separately by computing the relation 2]for a range of scaling 
factors k. From this follows the relationship for k in terms 
of a fixed desired magnitude M. 

This process provides a set of extinction-corrected SN 
la spectra that represents the most homogeneous spectral 
set we are able to generate. We now test these spectra for 
diversity in the sample. 



3 MEASURING SN lA DIVERSITY 

For each object we compile a single spectrum based on the 
average of all available spectra for that object between —2 
and 2 days from maximum. We use logarithmic bin spac- 
ing to offset the Doppler broadening of linewidths at higher 
wavelengths. The eight average spectra are binned to con- 
stant resolution, R = A/AA, with bin size proportional to 
the bin centre, 

(.) 

^ A refinement of this technique is to fit to both Mb and My 
simultaneously, which creates a wider, more shallow dependence 
of variability on wavelength band, further reducing the intensity 
relative to the spectral variations. This is achieved by minimising 
a best-fit function for fixed magnitudes Be and Vc, 

e{k) = {B-Bcf + {V -Vcf (1) 
= {a\ogk + h-Bcf + {c\ogk + d-Vcf . (2) 

Differentiating with respect to log k and setting the result equal 
to zero yields 

. = exp{ "(^-J+;f -'^) }. (3) 

The values of the coefficients are found as in the single-band case. 
Equation |3 allows us to determine a priori the scaling factor 
required to fix a spectrum near a chosen magnitude in both B 
and V. 
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Table 2. Catalogue of photometric parameters for SNe la. 



SN 


Ami5(B) 


-^max 




E{B - V) 


Mo 


References^ 


1981B 


1.10(07) 


12.00(02) 




0.11(03) 


30.10(30) 


Br83 


1989B 


1.31(07) 


12.34(05) 


12.02(05) 


0.34(04) 


29.73(04) 


We94 


1992A 


1.47(05) 


12.60(20) 


12.61(05) 


0.00(02) 


31.59(05) 


Ph99; Le93 


1994D 


1.32(05) 


11.84(05) 


11.92(05) 


0.00(02) 


31.10(07) 


Pa96; Ph99 


1994S 


1.10(10) 


14.85(04) 


14.84(06) 


0.00(03) 


34.36(03) 


Kn03; Ph99 


1996X 


1.31(08) 


13.24(02) 


13.21(01) 


0.01(02) 


32.19(07) 


Ph99; SalOl 


1998aq 


1.14(10) 


12.39(05) 


12.57(05) 


0.014(05) 


31.72(14) 


SahOl; Br03 


1999ee 


0.94(06) 


14.93(02) 


14.68(15) 


0.30(04) 


33.55(23) 


Ha02; St02 



^References.— Br83: Branch et al. 1983; Br03: Branch et al. 2003; Ha02: Hamuy et al. 2002; 
Kn03: Knop et al. 2003; Le93: Leibundgut et al. 1993; Pa96: Patat et al. 1996; Ph99: Phillips et 
al. 1999; SahOl: Saha et al. 2001; SalOl: Salvo et al. 2001; St02: Stritzinger et al. 2002; We94: 
Wells et al. 1994 



Our choice of /3 = 0.005 corresponds to approximately 20A 
bins near 4000 A and 35 A bins near 7000 A. This choice of 
binning parameter is narrow enough to appropriately sam- 
ple all prominent spectral features, but wide enough to be 
insensitive to noise, allowing for effective synthesis of the 
spectra. We tested a variety of bin widths from [3 — 0.001 
to /3 = 0.02 and the results are not sensitive to the width 
chosen. 

Using the averaged spectra from the eight supernovae 
in our sample we calculate the overall mean spectrum and 
determine the root-mean-square (RMS) deviation from the 
mean in each bin to quantify the diversity in that region. The 
result is plotted in Figures Q and |21 where the lower panel 
shows the overall average spectrum and the solid line in the 
middle panel shows the RMS dispersion about the mean. 
The shaded regions in this figure represent the uncertainty 
in the RMS dispersion, which is estimated by bootstrap re- 
sampling (Efron 1982). 

We wish to isolate the intrinsic diversity of SNe la from 
observational scatter. As the spectra used in the sample 
range from epoch -2 to 2 days we expect some diversity 
within the spectra for a single object. Moreover, differences 
between individual observations introduce scatter that is in- 
dependent of intrinsic diversity in the SN la population. We 
seek to quantify this variance and remove it from the to- 
tal, leaving only the intrinsic RMS. We have already made 
a single spectrum for each object by taking the average of 
all its available spectra. For each object, we account for the 
spurious dispersion by calculating the RMS variance of the 
available spectra for each object about the mean of that 
object. This quantity — the intra-object variation — when 
averaged over all the objects, gives an estimate of the dis- 
persion that is not due to intrinsic differences between the 
supernovae. This quantity is also used as an estimate of the 
non-intrinsic dispersion for an object when only one spec- 
trum is available. This dispersion is the dashed line in the 
middle panel of Figures Q] and It is small compared to 
the dispersion between objects. Since the contribution from 
observational error is expected to only increase the scatter 
in the results we subtract (in quadrature) the observational 
scatter component from the overall RMS deviation. The re- 
maining RMS, ideally all due to intrinsic diversity, is shown 
in the upper panel of Figures Q] and |21 

This method allows us to separate spectral diversity 
from scatter in absolute magnitudes. It will be sensitive to 



diversity in spectral features as well as colour variaions — 
large scale spectral shape differenes. The band to which we 
scale the spectrum will have the least diversity because we 
are effectively anchoring the dispersion to be lowest at this 
point. So when we are scaling to a fixed Mb the scatter is 
smallest between 4000 A and 5000 A, whereas when scaling 
to Mv the scatter is smallest between 5000 A and 6OOOA. 
The change in the overall shape of the dispersion spectrum is 
slight, and the significant spectral features remain unaltered. 
In Figure El we compare the final results for scaling to the V 
band (upper panel) and B band (lower panel). 

In Fig ure m we show a typical SN la spe ctrum at maxi- 
mum light ('Nu gent. Kim PerlmutteI^2002^ . The spectrum 
peaks at 4000 A. Our analysis shows that blueward of this 
peak the UV edge of the spectrum shows more variety than 
the rest of the spectrum. Much of this is likely to be intrin- 
sic SN la diversity, but it is also the region most susceptible 
to extinction and to observational flux errors. Rejecting the 
blue end of SN la spectra will reduce the scatter in SN la 
magnitudes that occurs both due to intrinsic SN variability 
and imperfect extinction correction. The Ray leigh- Jeans tail 
to the right of this peak, is a very stable region of the spec- 
trum until the Si II feature around 6 150 A produces a peak 
in variability. Beyond the Si II feature the diversity also ap- 
pears to increase. However, this may be another indication 
of imperfect extinction correction. 

We expect the dominant source of observational error 
that cannot be accounted for by measuring intra-object dis- 
persion to be due to extinction correction. By sampling the 
averaged object spectra with replacement many times we 
estimate the error in the extinction correction without mak- 
ing any assumptions about the error in colour excess values. 
Because we anchor the spectra to the B or V magnitudes 
of the supernovae (even when we are anchoring to the B or 
V magnitude of each supernova individually) , any faulty ex- 
tinction correction would show up as a bowl-shaped curve 
on our RMS vs wavelength plots (affecting high and low 
wavelengths equally - even though extinction preferentially 
affects short wavelengths) . The curves in Figure El may be 
the most illustrative for this point, since we have scaled the 
spectra to a fixed magnitude in B and V (upper and lower), 
effectively minimising the scatter in these regions. An error 
in extinction correction shows up as a gradual increase in 
the scatter away from these regions. 
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Figure 1. To calculate the diversity in SN la spectra as a function of wavelength we first calculate the average spectrum (bottom panel) 
and the RMS deviation from this average (middle panel, solid line). The uncertainty in this RMS is approximated by bootstrap resampling 
and is shown by the grey shading. The dashed line in the middle panel is the intra-object dispersion, and represents our estimate of the 
contribution from observational uncertainty. The upper panel shows our final result, which is the RMS minus this estimated observational 
dispersion. In this plot the spectra have been scaled to a fixed magnitude in the B band. This artificially lowers the scatter in the B-band 
region, centred around 4200 A. Large-scale shape variations in the spectra, would show up as a slow increase in the scatter away from 
this region (as we see here). Errors in extinction correction would also have this effect. 



In this analysis we are particularly concerned with (and 
sensitive to) small scale variation such as the variabiliy of a 
particular spectral feature. Prominent spectral features cor- 
respond to regions of high scatter, most noticably around 
the Si II line at 6150A. Ideally we would avoid these regions 
when calculating the magnitudes of SNe la for use as stan- 
dard candles. However, a large variation in a spectral line 
does not necessarily translate into a large variation in broad- 
band magnitude if the regions around the line are bright and 
stable, in which case the variability of these lines will find 
the greatest utiUty as a tracer of different subsets of the SN 
la population. 

In any real experiment we have a limited filter set, so 
cannot choose to observe the ideal region of every SN spec- 
trum for the continuum of redshifts we need to cover. How- 
ever, we can weight the error bars of our observations to 
favour those objects whose redshifts fortuitously allow one 
of our filters to cover the most stable region, and we can 
design filters narrow enough, and with sufficient overlap, to 
maximise the chance that a filt er will cover only the most 
stable region (|Davis et al.ll2nn6 h. 

Our analysis does not take into account any spectral 
differences between supernovae with different Amis values. 
When enough data is available it will be straightforward to 



apply this analysis to identify spectral features associated 
with wide or narrow light curves. 



4 MEASURING PECULIARITY 

When using SNe la for cosmology we need to select the 
observed supernovae with the most homogeneous character- 
istics for use as standardisable candles. Most well observed 
supernovae exhibit some peculiar features, and their full util- 
ity in cosmology requires the characterisation of these dif- 
ferences. A significant step toward this is the creation of an 
objective measure of peculiarity that can act as an objective 
criterion for including an object in a sample population, and 
also identify subsets of SNe la with very similar properties. 
Here we outline two different methods for quantifying pecu- 
liarity and demonstrate them on the small data set currently 
available to us. These methods will benefit greatly from a 
much larger \ow-z spectroscopic sample, such as is now be- 
ing collected by a number of groups, including the Nearby 
Supernova Factory, the Canada France Hawaii Telescope 
Legacy Survey (CFHTLS), the European Sup ernova Collab- 
oration and the Carnegie Supernova Project (|Hamuv et al.l 
I2OO5). In addition, applying these methods to the high-z 
samples being collected by the ESSENCE collaboration, the 
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Figure 2. As for Figure Q though here the spectra are scaled to have a the same magnitude in the V band (centred around 5200A), 
rather than the B band, as can be seen from the wavelength range of the region of lowest diversity. The same small-scale spectral features 
remain prominent, demonstrating the independence of the variance of these regions on the choice of anchor filter. For example, the 
variable strength (and position) of the Si II line observed at around 6150A. 



CFHTLS and the SN Cosmology project will allow for tests 
of evolution in the demographics of the SN la population. 

The first method uses a 2;-test to determine the degree 
to which a SN la deviates from the norm, and is particularly 
good at finding spectra with variations in regions which are 
usually stable, or spectra with strange continuum profiles. 
This means it is also good at picking out dubious extinction 
corrections. The second method uses principle component 
analysis and is particularly good at finding subsets of the 
population with similar characteristics. For example, it can 
pick out SNe la with high or low ejecta velocities by empha- 
sising correlated differences between spectral features. 



4.1 Peculiarity Parameter 

We use our calculation of intrinsic diversity to give a measure 
of the range within which supernovae can be considered nor- 
mal at each wavelength. If an object deviates significantly 
from the average in a region of the spectrum which is very 
stable, such as the range from 4000 A to 6000 A, we might 
consider it 'peculiar'. This can be quantified by assigning 
a z-score in each wavelength bin and then averaging the z- 
scores over some region, e.g. the entire spectrum or a band 
such as B or V. The z-score at wavelength A is, 

zx = (6) 
a 



where /, / and a are the flux of the object, the flux of the av- 
erage spectrum and the dispersion of the average spectrum, 
all at the wavelength A. We define the 'peculiarity parame- 
ter', P, to be the absolute average of these terms over the n 
bins in the desired wavelength range. 




A 



Subsets of this range, weighted by filter value, can be used 
to define the broadband peculiarity parameters, e.g. for 5, 

Pb ^ — y^^Ai^Ai, 

Z^A^A ^ 

where t\ is the throughput of the B band filter at wavelength 
A. 

Figure |5] shows a histogram of P for the objects in our 
sample calculated according to this method using the wave- 
length range 3700 to 7100A, based on the object averages 
for which the spectra were anchored in the B band. For this 
test, the average and RMS spectra are recalculated including 
observations from SN 1991bg and SN 1991T, previously clas- 
sified as 'peculiar'. The average spectrum changes negligibly 
when these objects are added. The peculiarity of 1991bg is 
made obvious by the test, from which we conclude that the 
spectrum of this object departs markedly from the mean in a 
non-uniform (across A) manner. However SN 1991T, which is 
known to be bolometrically over luminous, shows little spec- 
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Figure 3. Intrinsic RMS (total RMS minus observational dispersion) for spectra normalised to the average absolute magnitude in B 
(upper) and V (lower). 



tral peculiarity when scaled in magnitude to match other 
events. 

To test the fidelity of this parameter, the histogram is 
compiled again using instead the spectra anchored in the V 
band. Figure shows the recompiled histogram. Some slight 
rearranging occurs — notably the objects seem less bunched 
— however the prominence of 199 Ibg, and the relative lack 
of prominence of 199 IT, are unchanged. We conclude that 
this parameter can distinguish between events that exhibit 
peculiar broadband and spectral features, such as 1991bg- 
like events, but does not distinguish objects like 199 IT that 
are peculiar at wavelengths that also show significant diver- 
sity in normal supernovae. 



4.2 Principle Component Analysis 

Principle component analysis (PCA) is a technique for iden- 
tifying common characteristics in data with many vari- 
ables (e.g. Francis and Wills 1999). PCA has been success- 
fully applied to astronomical spectra i n the past to identif y 
different types of active galactic nuclei (jPrancis et al.l[l993 ). 
Here we demonstrate the concept on type la supernova spec- 
tra. This technique may be used in the future, with a larger 



sample set, as another method to identify subclasses of type 
la supernovae with similar properties. This is a major aim 
of several ongoing low-redshift supernova surveys and is im- 
portant for the control of systematic errors in planned (and 
ongoing) high-redshift surveys that aim to measure the ex- 
pansion histo ry of our universe. W e use the PCA technique 
developed by (Prancis et al.l (|l992l) including their analysis 
code with some minor variations.^ 

We perform PCA on the corrected sample of normal 
supernovae as outlined in Table Q as well as extinction- 
corrected spectra for the peculiar SNe 1991T and 1991bg. 
The result is shown in Fig. Q where we have plotted 
coefficients of the first and second principle components, cal- 
culated for each supernova. This test successfully identifies 
these two supernovae as peculiar. Moreover, this method 
is able to distinguish 1991bg-like SNe from 1991T-like SNe 
because they each appear peculiar in a different principle 
component. This method therefore shows great promise for 
selecting sub-classes of SNe la. Future work will perform this 



^ Many thanks to Paul Francis for supplying us with this code 
and running the first analysis. 
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This work 




4000 5000 6000 7000 

Wavelength (Angstroms) 

Figure 4. The industry st andard for typical SN la spectr a 
are the templates provided in lNugent. Kim PerlmutteJ (|20o3)- 
These are updated periodically, and the update of April 2005 can 
be found on supernova.lbl.gov/~nugent/nugent_templates.html. 
Here we overplot our average spectrum on the Nugent template 
for a SN la at maximum light. The agreement is good. Our 
analysis shows that the UV edge of the spectrum, below the 
peak around 4100A, shows an increased diversity among SNe la. 
(Whether this continues below the limit of our useful data, at 
3700A, is an open question.) This also approximately corresponds 
to where our average spectrum starts to deviate significantly from 
the Nugent template - corresponding to the greater uncertainty 
about the typical spectrum of this region. We have shown that the 
Rayleigh- Jeans tail to the right of the peak is a very stable region 
of the spectrum, until we reach the Si II feature at 6150A. This 
is therefore the best region of the spectrum to use as a standard 
candle for cosmological parameter estimation. 




Figure 5. Histogram of peculiarity parameter, P. Spectra have 
been scaled to the average absolute magnitude in the B band. The 
average spectrum from which deviation is measured is composed 
of average spectra from the objects listed excluding 1991bg and 
1991T. 
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Figure 6. Histogram of peculiarity parameter, P, with the av- 
erage spectrum additionally excluding 1999ee. Spectra have been 
scaled to the average absolute magnitude in the B band. 
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analysis on multi-epoch spectra, in which sub-class trends 
should be even more prominent. 



o 

^ -2 
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Figure 7. For each SN spectrum we plot the weights of the first 
and second components of a principle component analysis. The 
underluminous event, SN1991bg, is identified as peculiar in the 
first principle component, and the overluminous SN 1991T is the 
most peculiar SN in the second principle component. This is par- 
ticularly interesting because the spectra have been scaled to the 
same overall fiux between 3700A and 7100A, so the peculiarity 
that is apparent here is a difference in the shape and features of 
the spectrum, rather than simply a shift in the luminosity. 



5 CONCLUSIONS 

We have used published spectral data to calculate the diver- 
sity of type la supernova spectra close to maximum light as 
a function of wavelength between 3700 A and 7100A. While 
results in the regions toward the UV and NIR are sensitive 
to systematics of our calibration methodology, the diver- 
sity is particularly stable between 4100A and 6OOOA in the 
rest frame of the supernova. This is therefore the ideal re- 
gion for calculating SN la magnitudes for use as standard 
candles. There is an increase in supernova diversity below 
4IOOA. As we move blueward, extinction correction also be- 
comes more uncertain and therefore the region below rest 
frame 4100A should be avoided when using SN la as stan- 
dardizable candles. Observations of this region remain use- 
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ful for const raining theore t ical m odels of SN la explosions, 
for example iNugent et alJ (|l995f) showed the U band is a 
powerful tool for discriminating between different processes 
occurring in SNe la. Redward of the stable region the Si II 
line, typically observed at 6150A, shows a large diversity. 
This feature has previously been used as a tracer of light- 
curve w idth ([ Nugent et^_ay^ and to identify peculiar 
SNe la (Branch et alJ[l993al lbr). Our analysis confirms that 
it is a variable feature of ordinary SNe la. Future work will 
extend this study to spectra over a wider epoch range. 

We also demonstrated two methods for determining the 
pecuUarity of type la supernovae, firstly by comparing their 
spectra to the average SN la spectrum while taking into ac- 
count the typical variability of different spectral regions, and 
secondly using principle component analysis. Both methods 
show promise for use as objective tests of peculiarity. PCA 
in particular may be a good way to objectively define sub- 
types of the SN la class. 

The methods demonstrated here will benefit greatly 
from more spectroscopic data, not only more objects, but 
also extended wavelength coverage into the UV and IR. Al- 
though hundreds of supernovae have now been observed very 
well photometrically, the sample that have spectra published 
near maximum remains small. New data soon to be available 
from a variety of searches will enable a more robust analysis 
based on the methods demonstrated here. Such an analysis 
would aim to understand the nature, rather than the extent, 
of the diversity in the population. Potential sources include 
variations in expansion velocities, chemical composition and 
ex-situ factors such as evolution and gravitational lensing. 
Relationships between separate supernova events could then 
be drawn by a study of correlations in spectral residuals. 
This subject will be tackled in future work. 
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